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Abstract 

We present a full implementation of a quantum key distribution system 
using energy-time entangled photon pairs and functioning with a 30 km 
standard telecom fiber quantum channel. Two bases of two orthogonal 
states are implemented and the setup is quite robust to environmental 
constraints such as temperature variation. Two different ways to manage 
chromatic dispersion in the quantum channel are discussed. 

Since the birth of quantum key distribution (QKD) QIEIO], a lot of research 
and discoveries have been made ^ |S1 |7j leading today to demonstrated long 
distance QKD (8j i^j E| and even to commercially available quantum cryp- 
tography devices. However, these systems rely on faint laser pulses containing 
about 0.1 photon per pulse to guarantee absolutely secure keys. Therefore, 
only a fraction of the pulses will lead to effective bit transmission. Another 
consequence is that more than one photon are present per pulse with non-zero 
probability. True single photon sources |12[I13) seem to be a promising solution, 
but existing devices are not yet usable for out-of-the-lab systems. The most se- 
rious alternatives are systems based on entangled photon pairs |14l 1151 1161 [T7j , 
but no real- world long distance QKD using this approach has been reported yet, 
mainly for two reasons. First, the setup itself is relatively complex. Second, the 
spectral characteristics of the photon pairs sources imply that chromatic and 
polarization dispersion effects are not negligible and increase the error rate. 
Consequently, one has to find solutions if the aim is to deploy systems over fibre 
telecom networks. 

The goal of the present article is to demonstrate solutions to the chromatic 
dispersion issue for an energy-time entanglement based QKD system using a 
standard fiber quantum channel. We improved a setup previously presented 
by our group jl7| by using dispersion compensation or spectral filtering and 
extended the transmission range. 

Our setup is based on a Franson arrangement |18) of interferometers. A 
parametric down conversion photons pair source is located between Alice and 
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Bob. They both have an unbalanced Mach-Zender interferometer with photon- 
counting detectors connected at all outputs. When considering a given photon 
pair, four different events can be detected by both Alice and Bob. First, the 
photons can both propagate through the short arms of the interferometers. 
Then, one can take the long arm at Alice while the other takes the short one at 
Bob. The opposite is also possible. Finally, both photon can propagate through 
the long arms. When the path differences of the interferometers are matched 
within a fraction of the coherence length of the down-converted photons, the 
short-short and the long-long processes are indistinguishable and thus yield two- 
photons interferences, provided that the coherence length of the pump photons 
is longer than the path difference. If one records events as a function of time 
difference between detections at Alice and at Bob, 3 peaks appear (figure^. The 
central one is constituted by the interfering short-short and long- long events. It 
can be distinguished from the others by placing a time window discriminator, 
and is used to isolate nonlocal quantum correlation between Alice's and Bob's 
detections. On figure H we see that to allow the window discriminator to take 
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Figure 1: Schematic diagram of a Franson-type interferometers-source arrange- 
ment and the corresponding event frequency plotted as a function of the dif- 
ference of detection time at Alice and at Bob. The resulting peaks can be 
approximated by Gaussian with FWHM Ar and are separated by a time AT 
corresponding to the path length difference. A window discriminator of width 
W is placed around the central peak. The two phases 0a and (ps modulate the 
correlations. 

the maximum of the central peak into account (to have an optimal detection 
rate) but at the same time avoid the side peaks (non-correlated side peaks 
detections introduce errors), the separation between the peaks AT must be 
significantly bigger than At. At is the RMS of all the temporal spreading 
contributions: photon's coherence time, electronic and detection jitter and, for 
dispersive medium between the source and the interferometers, the spreading of 
the wave packet due to chromatic dispersion. This condition reads: 

At < AT < tP"'"^ 

where tP"™P is the coherence time of the pump laser. 

Dispersive media introduce chromatic dispersion oi D = per unit length. For 
standard telecom fibers and light around 1550 nm this value is D = 17 ps nm~^km 
The gaussian uncertainty in time due to chromatic dispersion only is then given 
by ATdisp — D AA L (where AA is the photon's spectral width and L the length 
of the fiber). Consequently, the inequality 
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has to be fulfilled to perform QKD with negligible errors due to chromatic 
dispersion. To match this condition for a given length of fiber, one can use large 
AT, reduce the dispersion D, or use photons with fine spectral width. We do 
not have to deal with similar problems caused by polarization mode dispersion, 
as information is not encoded in polarization. 

Theoretically, AT could be as large as needed by increasing the path length 
difference of the interferometers. For example, a difference larger than 3 m is 
needed for 30 km of standard fiber with our source. However, interferometers 
of this size are hard to stabilize. Moreover fiber interferometers suffer from a 
chromatic dispersion difference between the 2 arms that reduces the visibility. 
Another point is that even if there is no overlap between the peaks, dispersion 
broaden them. Consequently, wider coincidence windows are required, increas- 
ing the error contribution of the detectors' noise. For theses reasons we did not 
implement this solution. 

In order to reduce D it is possible to use dispersion shifted fiber for the quan- 
tum channel as in ^2] • This solution is not relevant because the resulting setup 
would not be deployable in a real-world telecom network, as there are only few 
installed lines using this kind of fibers. 

Another possibility is to use dispersion compensation ^Sj, with the only draw- 
back of the added loss on the quantum channel (which lowers the signal over 
noise ratio). This solution is the first that we investigate in this article. 
Alternatively, we reduce the AA of the 1550 nm photons (for 30 km the above 
necessary condition implies AA < 5.5 nm) by inserting a bandpass filter on Al- 
ice's side. As the central wavelength and the spectral width of the 2 down 
converted photons are related by energy conservation, filtering on one side re- 
duces the spectral width of the twin photons detected on the other side. Bob's 
detectors are gated upon Alice's detections thus, even if the key rate is reduced, 
the signal over noise ratio remains constant. 

To implement the full BB84 protocol, two different measurement bases are 
needed. This can be done by using 2 interferometers or a fast switch inside 
a single interferometer. Here we use a birefringent interferometer where the 
phase applied on the photons depends on the polarization. Thanks to this, 
we implement the full energy-time entangled BB84 protocol |21 IB) with only 
2 interferometers instead of 4. More details can be found in 17 . Pairs of 
810nm~1550nm photons are produced from a type I configuration of KNbOa 
crystal pumped by a 532 nm continuous laser. The 810 nm photons are collected 
in a single mode fiber and sent to the Alice's interferometer which is made of 
bulk optical components. It has 4 outputs, each of them consisting of a single 
mode fiber coupling system followed by a passively quenched silicium photon 
counter (EG&G). The 1550nm photons are collected in a single mode fiber and 
launched into the quantum channel which is a spool of several kilometers of stan- 
dard fiber. It can be optionally followed by our dispersion compensation device, 
i.e. a spool of negative dispersion fiber. The Bob's interferometers is connected 
at the other side of the quantum channel, together with a polarization beam- 
splitter and Faraday mirrors system that are part of the BB84 implementation. 
Two gated-mode InGaAs detectors (idQuantique) are connected at the inter- 
ferometer's outputs. The path difference of both interferometers correspond to 
Im optical length, and consequently AT =3.3 ns. 

An electronic system is used to trig Bob's detectors whenever a photon is de- 
tected by any of Alice's detectors. This electronics can also be used to char- 
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Figure 2: Schematic diagram of the experimental setup (PBS: polarizing beam 
splitter; LPF: low-pass pump filters; IF: optional 2nm bandpass filter; DCF: 
optional dispersion compensating fiber; FM: Faraday mirrors; ABE: adjustable 
birefringent element used as fibre A/4 plate) 



acterize the system in real-time: the information about which of the Alice's 
detectors registered a count can optionally be coded and sent to Bob's through 
a synchronized classical channel (consisting of another spool of standard fiber) 
using a 1550 nm laser. Upon detection at Bob's side, this information can be 
used for sifting and to verify the bits. This electronics thus enables to have an 
immediate information about the sifted raw key creation rate and the error rate. 
The most relevant experimental parameters are the following: spectral FWHM 
of the down converted photons: 6.9 nm at 1550 nm, 2nm at 810 nm; probability 
of having a photon coupled into the quantum channel whenever the 810 nm sili- 
cium detector fires: 0.5; losses of Bob's apparatus: 5.4 dB; quantum efficiency 
of Bob's detectors: 10%; detection gate width W — 1.1ns; false counts per gate 
due to the average noise of both detectors at Bob: ~ 1.13 x 10~^; false counts 
per gate due to noise on Ahce's side: negligible; 

The whole setup requires very careful tuning. The most difficult part is the 
alignement of Alice's bulk interferometer. In particular, we have to pay atten- 
tion that the probabilities for the short-short and long-long events are equals for 
all combination of detections at Alice and at Bob [2111 ■ Bob's apparatus must 
also be aligned to ensure that the polarization transformations in both arms of 
the interferometer are identical. 

We used a quantum channel consisting of a spool of 31 km of standard fiber, 
inducing losses of 8.3 dB. Two different configurations were implemented to limit 
the effect of chromatic dispersion: 

In the first one, the dispersion compensating spool (OFS) was connected in 
between the quantum channel and Bob's apparatus. This device compensate 
a dispersion of Dcomp — 506psnm~^, which corresponds to about 30 km of 
standard fiber. It induces losses of 2.9 dB, and a delay corresponding to 4 km 
of standard fiber. A 31-1-4 km spool of standard fiber was consequently used 
for the classical channel. In this case we implemented the full two bases BB84 
protocol using 4 detectors on Alice's side. The total count rate at Alice was 
about 79 kHz. The visibility of the interferences was about 89%. 
In the second configuration, an interferometric bandpass filter was placed in the 
source apparatus just before the collection lenses of the 810 nm output port (see 
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IF on fig. I2J). We used a filter of 2nm FWHM centered at 814 nm (as tliis fias 
been determined to be optimal for production/collection) and consequently the 
central wavelength of the co-detected photons was 1536 nm. As the Gaussian 

2 nm filter acts on an approximately 2 nm wide Gaussian spectrum, we obtained 
a width of about 1.45nm FWHM at 814nm, and 5.2nm at 1536nm. The 
number of 810 nm photons detected was reduced by a factor of about 3. To 
facilitate the alignements and the measurements in this second configuration, 
we increased the coincidences count rate by using only one of the two BB84 
bases. The total count rate at Alice was about 36 kHz and the visibility of the 
interference was about 92% in this case. 

Figure shows the resulting temporal distributions of the events obtained 
for 4 setups. Using the first plot, we obtained the total electronic/detection 
jitter of value 0.7 ns, because it is directly equal to the FWHM of the (central) 
peak; we also use this plot to verify that AT = 3.3 ns. The second plot clearly 
shows that a dispersion reduction method is necessary. Indeed, in the absence 
of any dispersion reduction, the two side peaks overlap to a large extend the 
detection window. We determined that in this case, the error rate due to the 
only contribution of uncorrelated events counted inside the detection windows is 
already about 10% of the total count rate. Moreover, the peaks are more than 

3 time larger than the detection windows, lowering the detection rate. The 
effect of the dispersion reduction is clearly visible on the third and fourth plots, 
corresponding to the two configurations described above. Note that the available 
filter of 2 nm does not completely remove the effect of chromatic dispersion. 
However, a narrower filter would decrease even more the total count rate. 

We achieved key distribution for our solutions. Because of phase instabili- 
ties in the interferometers the duration of a key exchange was limited to about 
40-50 minutes, but this issue could be addressed by using actively stabilized 
interferometers [2Tj . During this period the quantum bit error rate (QBER), 
the ratio of the error rate over the total rate after sifting, was about 10% on 
average. Tabled summarizes the performance obtained in terms of the different 
sources of errors. The total QBER is the sum of several contributions: i) The 



Table 1: QKD performances for both implementations 



Configuration 


Sifted 


Opt. 


Acc. 


Detect. 


Disp. 


Tot. 




key rate 


QBER 


QBER 


QBER 


QBER 


QBER 


compensation 


23 Hz 


5.5% 


1% 


4% 


^0% 


10.5% 


filtering 


12 Hz 


4% 


1% 


1.7% 


0.5% 


7.2% 



optical error rate due to the imperfect contrast of the interferences, ii) The 
accidental coincidences due to the non-zero probability of creating two pairs 
during the detection gate time-width: the photons from two different pairs are 
not entangled and thus have a 50% probability of producing incorrect bits, iii) 
The detectors' noise, which is independent of distance. As the detection count 
decrease with the losses in the line this detector QBER contribution increases 
with the distance and is generally the main contribution for large distances, iv) 
The dispersion which makes some non-correlated events to be registered inside 
the detection gate, as explained above. To estimate this QBER contribution we 
used the measurement presented in figure 13 Using the fitted Ar of the peaks. 
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Figure 3: Distribution of coincidence detection times for different setup. The 
squares show the experimental data. The plain curves show the individual 
gaussian peak; the heavy curve shows the sum of the 3 peak's contributions 
which is fitted to the experimental data to extract the FWHM At of the peaks; 
the dashed lines show the position and width of the detection gate. 
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it is possible to numerically integrate the 3 gaussian curves inside the discrimi- 
nator window to obtain the part of detections which comes from the side peaks. 
This is the error rate arising from the chromatic dispersion effect. Note that, for 
given AT and W, this integral value grows rapidly with At when this variable 
reaches about one third of AT. The chromatic dispersion QBER is thus very 
sensitive to spectral width and fiber length. 

The detector QBER is more important in the compensation configuration mainly 
because the compensating device add losses on the quantum channel, reducing 
the signal over noise ratio by a factor 2. However, one of our Bob's detectors was 
about 5 times less noisy than the other. With two such detectors, the detection 
QBER for the compensating solution would drop below 1%. 
The numerical integration of the events peaks was also used to calculate the 
fraction of the incoming photons which are accounted for bits of key. If all the 
photons that are part of the short-short /long-long events were inside the detec- 
tion gate this number would be 0.5. In our case it is 0.38 for the compensation 
solution and 0.22 for the filtering solution as the events are more widely dis- 
tributed in the second case (see figure EJ. These values are taken into account 
as a small added loss for the filtering solution with respect to the compensating 
one. Table|21summarizes the different factors leading to the registered sifted key 



Table 2: Factors leading to sifted key rates 



Configuration 


Singles 
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Tl 


Tb 


Tc 


Vd 




9s 


Sifted 




rate 




dB 


dB 


dB 








key rate 


compensation 


79 kHz 


0.5 


8.3 


5.4 


2.9 


0.1 


0.38 


0.7 


23 Hz 


filtering 


36 kHz 


0.5 


8.3 


5.4 





0.1 


0.22 


0.7 


12 Hz 



rates, starting from the singles rate at Alice. These factors are: ^: probability of 
having a photon coupled into the quantum channel whenever the 810 nm silicium 
detector fires; T^: loss of the quantum channel; Tg: loss of Bob's apparatus; 
Tq: loss of the dispersion compensating fiber spool; rjd- quantum efficiency of 
Bob's detectors; rjg: fraction of the incoming photons which are accounted for 
bits; Qs'. proportion of the bits that remain after sifting. This number should 
be 0.5 for a perfectly balanced two bases system. In the compensation case the 
value 0.7 is explained by two reasons. First, the bases choice at Alice is biased 
by the differences between the four bulk-to-fiber coupling and detectors efficien- 
cies. Secondly, the passive choice of bases at Bob depends on the polarization of 
the incoming photons. These photons are only partially depolarized by 30 km 
of fiber, and the remaining polarization fiuctuate inside the quantum channel 
during the key exchange. In the filtering case, as we use only one base at Alice, 
we could achieve a value of 1 by tuning the passive choice at Bob if the incom- 
ing photons where perfectly polarized, but this is not the case. Moreover, as in 
the compensation case, the polarization fluctuates inside the quantum channel. 
Note that a Qs 0.5 does not impair on the security of the scheme |22j. 

From theses results, we see that the choice of the best suitable dispersion 
reduction method is a matter of tradeoff between QBER and key rate values, 
and is different for each particular setup. In our case, numerical estimations 
show that an optimized compensation solution is better in term of key rate 
beyond 15 km, for a given QBER. However, the amount of negative dispersion 
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introduced must be calculated specifically for a given length of the quantum 
channel fiber. The main practical advantage of the filtering solution is thus 
that the system can be uniquely designed to be usable over a wide range of 
distances. The resulting key rate decrease can be compensated by pumping 
the source with a more powerful laser, or by using a more efficient photon pair 
source such as periodically poled lithium niobate waveguide |23j. When this is 
possible, filtering is more useful for real application. The only problem resulting 
from a configuration using a heavily pumped source filtered at Alice is the 
increase of accidental uncorrelated coincidence counts. This issue can be solved 
by filtering the photons at Bob with a corresponding filter. As both wavelengths 
are correlated, the only drawback is the non unity peak transmission of the filter. 
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Figure 4: Maximal spectral width that keeps chromatic and polarization dis- 
persion induced QBER below 1% as function of the quantum channel length. 
Plain curve: our energy-time entanglement implementation, with AT = 3.3 ns; 
dotted curve: possible optimal polarization entanglement implementation with 
PMD values of 0.1pskm~2 

The spectral width of the photons is also a limiting factor for polarization 
entanglement based QKD. QBER increasing in the case of polarization QKD 
is due to polarization mode dispersion (PMD) that depolarizes low coherent 
photons. Figure 01 shows the calculated maximal spectral width of 1550 nm 
photons that keeps the chromatic/polarization dispersion induced QBER below 
1%. The curve for energy-time entangled photons is calculated using our partic- 
ular experimental parameters, while the one for polarization is calculated using 
QBER = 0.5 X (1 — DOP) where DOP is the average degree of polarization 
computed numerically from formulas developed in j24| , using the standard PMD 
value of 0.1pskm~2. We see that for distances up to about 100 km, the energy- 
time solution is more robust to large spectral width. For longer ranges both 
solutions become similar with a slight advantage for the one using polarization, 
because of the square root dependance of the PMD with distance. However, 
these curves strictly apply to dispersion QBER and do not take into account 
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technical difficulties related to the necessary active polarization state control. 

In this article, we presented two practical means of dealing with chromatic 
dispersion induced problems in long-range QKD using entangled photon-pairs: 
dispersion compensation and reduction of the photons' spectral width at Alice. 
Optimal parameters where investigated and these solution were demonstrated by 
implementing a partial and complete BB84-like protocol with a setup featuring 
characteristic that could lead to a real- world telecom applications. A secret key 
was distributed over 30 km of standard fiber at a sifted bit rate of more than 
20 Hz and with an average QBER below the 11% limit for absolute security as 
stated in |2S]. Higher key rate is possible using more efficient sources. However, 
for a practical implementation, actively stabilized interferometers are needed. 

The authors would like to thank Claudio Barreiro and Jean-Daniel Gautier 
for technical support. Financial support by the Swiss NCCR Quantum Photon- 
ics is acknowledged. 
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